CU-TP/11-03 



Quark/lepton mass and mixing in S3 invariant model 
and CP-violation of neutrino 

T. Teshima and Y. OkumuraQ 

Department of Natural Science and Mathematics, 
Chubu University, Kasugai 487-8501, Japan 

Abstract 

Weak bases of flavors (u, c), (d,s), (e,u), (ve^fi) are assumed as the S3 doublet and t, b, r, 
v T are the S3 singlet and further there are assumed 53 doublet Higgs {B.\,H.2) and S3 singlet 
Higgs Hs- We suggest an S3 invariant model in which the Yukawa interactions constructed 
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sector, neutrino masses are assumed to be constructed through the see-saw mechanism from the 
Major ana mass. The tri-bimaximaklike character of neutrino mixing Vmns can be explained 
dynamically without any other symmetry restrictions. It is predicted that a quasi-degenerate 
mass spectroscopy of neutrino is favorable, and values of | Vmns 1 13 > CP violation invariant measure 
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I. INTRODUCTION 



In present elementary particle physics, the problem of the origin of quark/lepton mass and 
mixing is the most interesting and challenging subject, for the exploring this problem leads 
to the finding of a clue of a new theory over the standard theory of the elementary particle 
physics. Although the quark mass hierarchy and quark mixing Vckm have been explained 
successfully by many authors but neutrino mass hierarchy and mixing Vmns having the 



large mixing character 1CH12|| is not explained sufficiently Especially, smallness of neutrino 



mass and tri-bimaximal-like mixing nature of neutrino are not explained satisfactorily in the 
same footing as discussion of quark mass and mixing. 

In these circumstances, many models in which flavors of quark and lepton are governed by 



the symmetry group 53 



13M22| , S4 or A4 23M25I] have been analyzed. Fundamental scenario 



to these models is that quark and lepton flavors and further Higgs fields are considered 
as to be governed by the discrete symmetry, S3, S 4 or A 4 group, and physical neutrinos 
are Majorana neutrino induced from see-saw mechanism through the mixing with right 
handed Majorana neutrino. Furthermore, almost models except ours 19[ have considered the 
additional constraint, Z2, S2 symmetry or fi — r symmetry, for explaining the tri-bimaximal- 
like mixing nature of neutrino. In contrast, we explained this tri-bimaximal-like mixing 
nature by a dynamical mechanism, in which neutrino mass matrix is induced from the 
mixing between Dirac mass and Majorana mass of neutrino, and the hierarchy between 
masses of S3 singlet and doublet. 

In our previous paper Jjj] , we used a standard Yukawa interactions modified in order to 
make mass matrices for quark and lepton Hermit. But this modification of standard Yukawa 
interaction violates the conservation of hyper charge Y, then in the present work, we use 
a standard Yukawa interaction without any modification. Further, CP phases, one Dirac 
and two Majorana phases in neutrino mixing, are analyzed and CP-violation measure J and 
effective Majorana mass | <m> | in neutrino-less double (3 decay (((3(3)ou) are estimated. 



II. S3 INVARIANT MODEL 



First, we explain our S3 invariant model. We assume that S3 symmetry governs the 
generations of quark and lepton (charged lepton and Dirac neutrino). Weak bases of flavors 
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(u,c)l,r, (d,s) L;R , (e,/i) LjjR , (v e ,v^ Lfi are assumed as the S3 doublets and t LtR , b LjR , t L;R , 
v tLR are the S s singlet. Further there are assumed £3 doublet Higgs (Hi, H 2 ) and S 3 singlet 
Higgs H s . If we represent the fields of quark and lepton as /, our S3 model are composed 
of the following fields, 



S 3 singlet : fe R , H s , 
S 3 doublet : fi R =(rf' R ,f? R ) T , H D =(H 1 ,H 2 f. 

and the fields f L ' R represent the following quarks and leptons, 



(1) 



\fs ) 



cL,R 
12 

f L,R 
JS 



UL,R 

cl,r 

\ t L,R J 



d>L,R 
SL,R 
\ h L,R ) 



v eL,R 

V t J -L,R 
\ V tL,R J 



I e L , R \ 

VL,R 

\ T L,R ) 



(2) 



\T 



In the SU(2)l gauge space, Higgs fields Hp, Hs are SU(2)l doublets; Hp = (Hjy, H' D j , 
H s = (H s h ,H%) T , respectively. For quark fields, Q\ = (u L ,d L ) T , Q 2 = (c L ,s L ) T , Qs = 
{tL,b L ) T are SU(2) L doublets, and df = d R , d R = s R , d§ = b R , uf = u R , u R = c R , u§ = t R 
are SU(2) L singlets, and for leptons, L[ = (veL^L) 1 \ L 2 = (v^ L , Hl) T ■, Lg = (u tL ,t l ) t are 



SU(2) L doublets, and Zf = e R , l R = fi R , l§ = t r , v r = v eR , v\ 
singlets. 

We start from the standard Yukawa interaction, 



VRi 



v R = v TR &re SU(2) L 



-C f D = J2 FtjkQiHjdZ + V u ijk QfeH*u R + T l ijk LfHjl R + T^ k LfeH*u R } + h.c, (3) 

i,j,k=l,2,S 

where V(- k are complex interaction strengths and e is the 2x2 antisymmetric tensor. We 
assume the following S3 invariant mass Lagrangian for quarks and leptons under the spon- 
taneous symmetry breaking of vacuum (Hs) = (0,Hg) T , (H D ) = (0,H° D ) T , 



-C% = T d s l f^f R Hl + r d Jj L D f R H° + T%[(f?f? + tif?)m + (fif R - f 2 L f 2 R )H°] 

+rMH° D f R + 7^H° D T f R ) + h.c, 

for down-type quark and charged lepton (4) 

-r 6 » = ryJif R Hf + ryj L D f R Hf + i%[Cf?f* + + (J[f R - Jff*)n?\ 

+T u MH° D *f R + J^H° D * T f R ) + h.c, 

for up-type quark and Dirac neutrino 
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where we used the fact that the 5*3 doublet can be made from the tensor product of f D and 



fo ll3Hl5| as 



in 



tff? + f£h R 

~cL fR _~cL fR 
Jl Jl J2 J2 

These mass Lagrangian (4) are almost similar to ones proposed in our previous paper 
which Hi coupled to Jiff and Jjjf? in C D d ' 1 and coupled to Jfff 1 and fgf? in C D U ' U 
are interchanged. In our previous model, the mass matrices produced from the S3 invariant 
mass Lagrangian are Hermit, but the Yukawa interaction for these mass Lagrangian dose 
not conserve the hyper charge Y. 

From the SU(2) L gauge freedom of fields H s , fs' R , fi' R and /jf'^, we can choose the 
phase of fg' R as 

phase of H s — phase of / 4 L + phase of f R = 0, 
i = S, 1, 2 

phase of f R = phase of f R = phase of f R . (5) 



Then we can get the mass matrices M^i for down-type quark and charged lepton, and M UtV 
for up-type quark and neutrino Dirac mass as follows; 



-C f 



fLM f f R + h.c, f = d,u,l,u, 



I d,l 1 d,l iSo 

1 n{ + H2 e 



M, 



d.i 



A/x 



d,l ifa \ 



\ d,l 

A/i 2 ' e 



/ u.u 1 u.u - 



d,l 
Ml 



-102 



V 



\f4 
\f4 



n e 

/if e i(t>2 



Mi' ~th' e ~ 



d,l 

Mo 



Vu,u 
3 e 

u,u - 

M 3 e 
Mo 



-i<j>2 



where we use the following parameterization, 

Mo = ( r sl-^sl)sS, Ml = (r^! 1) n,22, M2 = ( r D2 1 -^2 l)ll,22, 

A/^2 = (r{, 2 |-ff?|)l2,21, M3 = ( r D3 1-^2 1)25,52, A/J 3 = (^3 1 -H"i|)l5,Sl, 



A 



(6) 



phase of H® — phase of H s , (f) 2 = phase of H° — phase of H s . 



Yukawa interaction strengths for the interaction including f D or H° D are considered to 
be very small compared to T s for the interaction including fg and H s , then /i 2 , /x 3 <C 



Hq. Thus the mass hierarchy of d-type and w-type quarks and charged leptons are realized 
in our model. 

For neutrino mass, we assume that there are very large Majorana masses constructed 
from the right-handed neutrinos, the existence of which are suggested in the 5*0(10) GUT, 
and from this Majorana mass we can get the very small neutrino masses through the see- 
saw mechanism [26j. We assume that the Majorana mass is constracted as S3 invariant 
containing only right handed neutrino v R = (v R 7 v R ) T , v R and has no Higgs field [15I ]. 



C 



_ IpM/, .R\T r -l 1 R 1 '^ V M(,.R\T r i-l,.R 4- h n 

> ' - 2 s ^ s ) ^ v s + 2 D ^ D > D+ ' 



(7) 



where C is a charge conjugation matrix. This Majorana mass term is expressed as 



M 



M 



-{v R ) T C- x M M v R + h.c, 
' Mi ^ 



M 



Mi 
M 



(8) 



using the next parametrization 



M = (rf) ss , Mi = (r^) 11)22 . 



If Mi <C Mo, as the case of Dirac neutrino mass, we can explain the tri-bimaximal-like 
mixing character of neutrino without any other symmetry restriction. We will discuss this 
problem in section 4, in detail. 



III. QUARK MASS AND MIXING 

In this section, we consider the quark mass and mixing in detail. We assumed the 
Yukawa interaction describing the masses for ci-type and u-type quark as Eq. (4) and the 
mass matrices for these Yukawa interactions were expressed in Eq. (6), 



-£{, = f L M f f R + h.c, f = d,u 

( 



f 



L,R 



( d L ,R ^ 

SL,R 

\ h ' R J 



Ma 



A/^e^ 1 /if 



A/^e^ 1 A/ife^ 1 \ 
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./ 



L,R 



CL,R 



A/^e"^ 1 



Vs e ~ 



/ 



Because mass matrices M d and M u are complex symmetric ones, these matrices are diago- 
nalized by the unitary matrix U and Vas (see Appendix) 

VjM d U d = diag[m d , m s , m b ], V d = U d *S d \ S d = diag[e^, e^, e**\, 
a d, fid. Id are arbitrary real constants, 
(dg, sS, 6£) T = C/j (d fl , s R , b R f, 8 % b?f = Vj (d L , s L , b L f, (9) 

VjM u U u = diagK, m c , m t ], K = 5 U = diag[e ia ", e^«, e*»], 

Pu, 7« are arbitrary real constants, 

K, eg, = Ul (u R , C R , t R f, «, C™, f™) T = K f («L, CL, U) T , (10) 

where (m^, m 2 s) ml) and (m^, m;?, ml) are eigenvalues of M d M d ^ and M n M u ^, respectively, 
and (cf™, s m ,b m ) and (u m ,c m ,t m ) are the mass eigen states for the weak basis (d,s,b) and 
(w,c, t), respectively. As a result, the CKM mixing matrix Vckm in the weak charged 
interaction is expressed as 



Vckm = VjV d = S u UjU d *Sj. 



(11) 



This Vckm matrix can be parametrized by three mixing angles 9 12 , #23, #13, and a CP- 
violating phase 5 after adjusting 6 arbitrary phases in S d , S u , where 1 of 6 phases can be 
settled as 0. Standard expression of this matrix is written as 

Vckm = 



C12C13 S12C13 Si3e 

-S12C23 - Cl 2 S23Sl3e l<5 C12C23 - Sl2S23Sl3e* 5 S23C13 

Sl2 s 23 ~~ Ci2C23Si3e l<5 — C12S23 — Si2C23Si3e i<5 C23C13 

where Sij— sin 9^, Cij = cos 5 is CP- violating phase. 



V 



(12) 



We calculate eigenvalues of quark masses and diagonalization matrices U analytically, 



under the assumption / u 1 ,/i 2 ,/i3 <C Ho and A 



\H$\ 



< 1. 
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<0> d-type quark: 

/if + (1 + A 2 )/if + 2(4 J (cos 2 2 + A 2 cos 2 00/if + A 2 sin 2 (0i - 2 )/if ) 
/if + (1 + A 2 )/if - (cos 2 2 + A 2 cos 2 00/if + A 2 sin 2 (0i - 2 )/if ) 



"id 



1/2 



1/2 



, (13) 



m b « /i , 
/ 



cos #d 
— sin 9 d e~ l9 z 



sin d e* e 3 
cos# rf 



M[T e 

M3 r —i<h. 



1 ^(-Acos^e^ 1 +sin^e^ 2 - ie 3) -^(Asin^e^ 1+ ^3 +cos^e^ 2 ) 1 i 

\ Mo Mo / 



(14) 



tan6> 



d _ /4sin(0i - 2 ) 

3 — d A 

[l\ COS 01 



tan 



A^/ /if cos 2 0i + /if sin 2 (0i - 02 ) 



COS 2 2 + A 2 COS 2 0l)/if 2 + A 2 Sin 2 (0i — 02)/if + /if COS 02 



<) w-type quark: 



m. 



.-,1/2 



/if + (1 + A 2 )/if + 2/iW(cos 2 02 + A 2 cos 2 0i)/if + A 2 sin 2 (0i - 2 )/i : 



4 2 ) 



.-,1/2 



/if + (1 + A 2 )//f - 2/i^ J (cos 2 02 + A 2 cos 2 00/if + A 2 sin 2 (0i - 2 )/if ) 



(15) 



m* « /i , 



Mo 

EjLpi<i>2 
Mo 



cos # u sin ^ u e jy 3 

— sin 9 u e~ td % cos # M 

^ ^(-Acos^e-^ 1 + sin^e-^-^s) — g|(A sin^ e -^ 1+ ^3 + cosfl M e^ 2 ) 1 J 

(16) 

/if sin(0i - 2 ) 



tan 6% = 



(l\ COS 0i 



. v Vf cos 2 0i + /if sin 2 (0i - 2 ) 
tant7 u = . 

V (COS 2 02 + A 2 COS 2 0l)/i" 2 + A 2 Sin 2 (0i - 2 )/i2 2 + K COS 02 
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The CKM matrix is also written analytically using the expressions Eqs. (14) and (16) as 



,t 



Vckm — VjVd — S u U u T Ud* Sd' 

• cos 9 U cos 9 d + sin 9 U sin 9 d e i{ ~~ e ^ +e ^ 



sin 6 U cos 9 d e %e % — cos 9 U sin 9 d e %e z 



11 -sin^e-^ 2+ie 3)4 

Mo 



, (Acos^e^ 1 - sin^e^ 2+ie 3)j| - (Acos^e 

cos 9 U sin 9de~ td z — sin 9 U cos 9 d e~ ld ^ 
cos U cos d + sin 9 U sin 9 d e ie 3~ ie 3 
(\sm9 d e i ^- i6 i ± cos^e^ 2 )4 - (A sin^e"^ 1 "^ + cos^e"^ 2 )4 

'Mo Mo 

(Acos^e^ 1 - sin 9^-^)4 - (Acos^e - ^ 1 - sin^e"^ 2 "^^ \ 

Mo Mo 

(Asin^e^ 1+i9 3 + C os0„e^ 2 )4 - (A sin 9 u e~ i<l>l+ie 3 + C osfl u e-^ 2 )4 

Mo Mo 

1 / 



(17) 



Next, we examine our model numerically. The present experimental data for quark masses 
and CKM matrix are given in the PDG 2008 (2t[ ; 

TO i fry 

— = 0.045 ± 0.025, — = 0.025 ± 0.008, m b = 4.20 ± 0.12GeV, 

m s m b 



m„ 



0.0019 ±0.0008, 



m t 



W< 



CKM | 



0.0074 ±0.0006, m t = 171.3±2.3GeV, 



0.9740 to 0.9744 0.2247 to 0.2267 0.0034 to 0.0038 
0.2246 to 0.2266 0.9731 to 0.9736 0.0404 to 0.0425 
0.0084 to 0.0090 0.0397 to 0.0417 0.9991 to 0.9992 



(18) 



, _ ±0.031 _ ±0.015 

vertex coordinate of unitarity triangle p = 0.135 , ri = 0.349 . 

J 6 y -0.016' ' -0.017 

±0.19 

invariant measure of CP violation J = (3.05 ) x 10~ 5 , 

v -0.20 7 

where vertex coordinate of unitarity triangle (p, fj) and Jarlskog invariant measure of CP 
violation J are defined as 



p = Re 



(Vckm) n (Vckm) 13 

_(VcKM)2l(Vc KM )23j 



7] = Im 



(VcKm)ii(VcKm)i3 
.(VcKm)2i(^CKm)23 



18') 



J — Im [(Vckm)i2(Vckm)23(Vckm)i3(K:km)22] • 
Using the computer simulation, we estimate the allowed region for values of 11 parameters 



(A*?)) A*2' A*3' A*0' A 4 !' A*2 ' A*3 ' ^ 



i> 2 ) so that the quark masses and Vckm computed 



from these parameters satisfy the experimental data Eq. (18). The estimated results are 
following values; 

t4 = 4.20 ± 0.12GeV, ^ = 0.0120 ± 0.0030, 4 = -0-0136 ± 0.0004, 

^| = ±(0.0282 ±0.0008), 

^ = 171.3 ±2.3GeV, ^ = 0.00369 ±0.00003, = -0.00378 ± 0.00003, (19) 

c| = ^(0.0127 ± 0.0007) (opposite sign for that of the ratio (4/f4), 
A = 0.207 ±0.004, 0x = -(74.9 ±0.8)°, 2 = (0.74 ± 0.31)°. 

For these values of parameters, three mixing angles and a CP-phase in standard expression 
Eq. (12) for Vqkm are estimated as follows 

012 = (13.1 ±0.1)°, 23 = (2.38 ±0.05)°, 13 = (0.207 ±0.009)°, 5 = (68.1 ± 3.8)°. (20) 

Finally, we consider what are the origin of the Cabibbo angle and CP violation phase in 
our model. |Vckm|i2 elements is expressed approximately in Eq. (17) as | cos 9 U sin 9 d e~ t0 3 — 
sin 9 U cos 9de~ ie % | . This is expressed by using the approximation 3 —0" ps 0i and 0^ ~ M , 
which are recognized from the expression for 0|, 3 in Eqs. (14) and (16) and the values Eq. 
(19), as follows 

| cos 9 U sin 9 c ie^ 1 — sin 9 U cos c ze i( ^ 1 1 = [(cos 9 U sin 0^ — sin 9 U cos 0^) 2 cos 2 0i ± 
(cos0 n sin0 rf + sin0 n cos0rf) 2 sin 2 0i] 1 ^ 2 ps | sin(0 d + U ) sin 0i| ps |Asin0i|. 

We can say that the origin of Cabibbo angle and CP violation phase are the ratio A = 
\H®/H$\ and the relative phase 0i between phase of if° and if<j. 



IV. LEPTON MASS AND MIXING 

In this section, we consider the charged lepton and neutrino mass hierarchy and neutrino 
mixing Vmns, which has the tri-bimaximal-like mixing character. We assume the Dirac mass 
terms Eq. (6) for mass of charged lepton (e, /i, r) and Dirac mass of (u e , u T ) and the 
Majorana mass term Eq. (8) for Majorana mass. For charged leptons (e, //, r), mass term 
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is represented in Eq. (6) as 



( °- - ^ / fi[ + /4e^ 2 A/4^ 1 A/ifje^ 1 ^ 



-C l D = l L M t l R + h.c, 
eL,R 

\ T L,R ) 



l L > R 



M, 



A/4^ 1 n\ - jjje^ 2 y^e** 2 



\ Xpt l 3 e^ 



l p i(f>2 



fi' 3 e 



This mass matrix M t is diagonalized by the formula similar to Eq. (9), 

VfMM = diag[m e , m^rrir], V x = UfSj, S t = diag[e^, e^, e*»], 
an, A? li are arbitrary real constants, 
(eS, r£) r = U! (e R , fi R , r R ) T , ( e ™, r?) T = (e L , ^, r L ) T , 



(21) 



where (m 2 e , m 2 ^, m 2 .) are eigenvalues of M/M/ t and (e m ,/i m ,r m ) are the mass eigen states. 

For the neutrinos (v e , v^, u T ), we assume the Dirac mass Eq. (6) and Majorana mass Eq. 
(8), 



— C v 

*~-D+M 



v L M v D v R - -{v R ) T C~ l M M v R + h.c, 



(22) 



where 



v = 



( \t\ + /^e"^ 2 X^e'^ 1 A/i^e"^ 1 N 



Vs e " 



K J 



f Mi ^ 



M M = 



M x 
M 



Using the relation, 



u L M v D u R = \(v L M v D v R + v L (M v D ) T v R ), 



where v is an anti-neutrino of v, and v R = Cu L , v L = —v R C~ l , we can rewrite the 
neutrino mass terms Eq. (22) to the 6x6 matrix as 



J -D+M 



1 

2 



(v L , L ) 







M 



•T 



D 




+ h.c 



(23) 
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Because M V D <C M M , this 6x6 mass matrix in Eq. (23) is diagonalized by the 3x3 matrices 

as 



J ~D+M 



'J)R V R\ 



-M" D M^Ml T 
M M 

-(M^M^M^Y 







M. 




The (1,1) block of the matrix in this equation, M^f = -M V D M^M V J is very small compared 
to the (2,2) block, M M - This result is the see-saw mechanism. We analyze the mass terms 
responsible for the mass matrix M^f = -M^M M l M^ T . 



'-D+M 



U L M M R + R M M< U 



Aft L 



(24) 



= -M» D M M X M V J. 



Mp is a complex symmetric matrix and M M is a real symmetric matrix, then M^f is a 
complex symmetric matrix. Then the M^f can be diagonalized by the unitary matrix U u 
and V v , as for charged leptons, 



V}MjfU v = diagK e , m u . m Vr ], V u = U* u , 



(25) 



j,m ,/m ,/m 
\ U e R-> u n R > u r RJ 



where (m£ e , ml , m 2 UT ) are eigenvalues of Mff Mff^ and u m = (v™, v™, v™) T are the mass 
eigen states. Here, it should be pointed out that the relation V v = U* in present neutrino 
case is caused from the fact that the Majorana mass term v L M^fv R = — v rT C' x M^fv R is 
consisted from only the v R field. Using the mass eigen states u m , we rewrite Eq. (24) to 



'-D+M 



-R 



R 



u mL diag[m Ue , m v , m Vi )v m + v m diag[ra„ e , m v , m VT \v 



,mL 



^X m diag[m, e , m Vfl , m UT } X m , X™ = » mL + 



-R 



V" 



(26) 



where x m is a Majorana neutrino. The mixing matrix of neutrinos corresponding to the 
Vckm in quark sector, Vmns is expressed as 



V f 



MNS 



(27) 



This Vmns is parametrized by three mixing angles and one Dirac CP-violation phase and 
two Majorana CP-violation phases, adjusting 3 arbitrary phase parameters in Si. Standard 
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expression of this matrix is written as 
/ 

Vmns = 



V 



C12C13 S12C13 

-S12C23 - Ci2S23Sl 3 e^ C12C23 - Sl2S2 3 Sl3e* 5 
■S12S23 ~ Ci2C23Si3e l<5 — C12S23 — s 12C23 s 13 e * 



id 



si 3 e 
S23C13 
C23C13 



\ 



/ 



Mi 



where Sj,- = sm#j,-, Cjj = cos 5 is CP- violating Dirac phase, 



(28) 



Pm — diag[l, e* 2 , e* 2 ], /3, 7 are CP-violating Majorana phases. 

Diagonalizing Mff , we calculate the eigenvalues of Majorana neutrino mass and neutrino 
mixing matrix. First, we calculate M^f = —M^M M 1 M^ ) T , neglecting A 2 and 2 because of 
the results that A ~ 0.21, <p 2 ~ obtained in the quark sector numerical analysis, 



M 



M 



( 



-M V D M^M V J 

Mi 



\(ftj±± 1 \^_\p-i<t>l \ ^3(^1+2^) 



'H Mi " r M / e 



MM. 

Mi 
Mi 



Mq ' 



" r M 



\ A v Mi M / Mi M 



Using the assumption Mi -C M and | /x^ 

( K 



Mi 

Mi ~ r M 

u u2 n" 2 

M-3 1 Mo 

Mi M 



^2 1 ~ I A*3 1 ^ I A*o 1 5 can ^ e expressed as 



/ 



M 



M 



1 

m7 



2A/x ! 1 '/x^e-^ 1 

(K - /4) 2 

^ A fa(fi + 2^) + ^g^) - ^) + 

A fate + 2M + 2*0) e~*\ 



M 



M 



1M0M3 
M 



/4 



M, 



(29) 



/ 



It is recognized that this mass matrix M^ 1 has a character that v^-v? mixing becomes 
maximal and f e -f M mixing can be large. This can be confirmed from an approximation to 
neglect the terms including M\ /Mo, and to set — ~ l/^l, Hi + = ^ ^ A 4 !? induced 
from the assumption /i^ —fi 2 ~ I A^l ^ Akb which is obtained in previous quark sector 
analysis. Thus the mass matrix Eq. (29) can be parametrized as 



M 



M 



1 

Mi 



5 2 



. A i/2 -j0i 
2^3 e 



2^3 e 
i/2 



V -^3(^3 -35)e"^ ^ 



/if 



(30) 
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This mass matrix can be diagonalized analytically by the unitary matrix U u and L/J as given 
in Eq. (25), 



U/M^U U « 



diag 



5 



d 2 



2& 



2M^^^ 2M^^ + ^ M x 



where a 
( 



y/2S' 
10 ^ 

l/y/2 l/y/2 
^0 -l/y/2 l/y/2 J 

where rj w -={ 

2y/2 



1 





r]e % ' 

1 



>t\ 



35 



/ 



\ 



\ 



-i<t>i 


1 


/ V 




tan 26 ? 






" \/2<5 


i- 

e 2 





\ 
















e 2 / 



cos# sin^e"^ 1 
sin 6e i(f>1 cos 9 
1 I 



S u , (31) 



From this approximate expression, it is recognized that v^-v T mixing becomes maximal and 

3 

v e -Vn mixing angle can be large, for example, if 5 ~ A//3, tan 2$ w — =. Furthermore, it 

\/2 

is recognized that |Vmns|i3 ~ V ~ A/2y / 2 is small but not 0. Thus, the tri-bimaximal-like 
mixing character of neutrino can be explained dynamically without any other symmetry 
than 5*3 symmetry. 

Next, we examine our model numerically. For charged lepton sector, experimental data 
of the masses of (e, /1, r) are given in PDG 2008 [27] as, 

(32) 



7TI 777 

— = 0.004836 ± 0.000001, — ^ 



0.05946 ± 0.00001, m T = 1776.84±0.17MeV. 



m. 



For the neutrino, the experimental data obtained from the neutrino oscillation is summarized 



in Ref. 



12|, as 



6.90 x 10" 5 eV 2 < Am 2 < 8.20 x 10 _5 eV 2 







0.27 < sin 2 # Q < 0.37, 



2.15 x 1(TW < \Ami t J < 2.90 x 10~ d eV 2 , 
0.33 < sin 2 # atm < 0.65, 
[VWshs 2 < 0.052, 

where these values are including 3a uncertainty. We assume that |Am 2 

2 m 2 . and ni „ 



(33) 



atm I 



\ Am ll\ 



|Am 32 |, Am Q = Am 21 , where Am^ = m~ v . 



9 2 3 mixing angle, 9 Q = 9 12 
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eV 

0.08 
0.06 
0.04 
0.02 
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- 
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- 




mv3 














































mv2 




Ivi 

















0.0 



0.2 



0.4 0.6 

ri/2=mvi/mv2 



0.8 



1.0 



FIG. 1: Values of function of ratio ryz = m vl /m U7 . This figure is obtained 



from the values of Am|, = Am^ and Am^ tm 



Am^ in Ref. 



12 ], assuming the normal 



hierarchy m vi < m U2 < m U3 . Upper two lines denote the m U3 values with 3a error, middle 
two lines the m U2 values and the lower lines the m Vl values with 3<r error. 



mixing angle. In our analysis, we assume a normal hierarchical (NM) mass spectrum 
m Vl < m„ 2 < m„ 3 for the neutrino masses. From the experimental data for mass squared 
differences Am^ tm , AttIq, we can obtain the numerical values of function 



of ratio ryi = m Vl /m h 



2 > ITlvi 
\/2 



ri/2in U2 , m U3 



shown 



vi ■ 



in Fig. 1. For the value near 1 of the parameter rx/2, the quasi-degenerate (QD) mass 



spectrum m v 



m U3 appears. From the experimental data for mixing angles, the 



magnitudes of elements of Vmns are restricted as 

/ 

I^MNsI = 



0.77 
0.10 
0.28 



0.88 0.46 
0.49 0.47 
0.61 0.34 



0.61 0.00 
0.78 0.57 
0.71 0.57 



0.26 
0.81 
0.81 



\ 



/ 



(34) 



Using Eqs. (6), (21), (22), (24), (25), (27) and the numerical result (19) for A, fa and 2 
determined in quark sector analysis, (32) for charged lepton mass and experimental data 
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0.0 u 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 u 

0.0 0.2 0.4 0.6 0.8 1.0 

ri/2=mvi/mv2 



FIG. 2: The ratio m„ 2 /mj, 3 as a function of the ratio ri/2 = m Ul /m U2 . Two lines are obtained 
from the experimental data (Ref. |12l]). Sequence of dots plotted in middle is obtained by 
satisfying the allowed region (34) of neutrino mixing | Vmns I i n the case of 0i = 0° , 02 = 0° , 
where 0i i2 axe included in mass matrices, Eq. (6). Sequence of dots plotted in right-hand 
side corresponds to the case of <f>\ = —74.9°, 02 = 0.74°, which are determined in quark 
sector analysis. 

(34) for neutrino mixing, we estimate the allowed values for parameters fj\ and Mf, 

A = 0.207 ±0.004, 0x = -(74.9 ±0.8)°, 2 = (0.74 ± 0.31)° 

(in-put data determined from quark sector analysis), 

fj} = 1776.84 ±0.17MeV, ^ = 0.0308 ± 0.0007, = -(0.0307 ± 0.0017), (35) 

A*o Mo 

= -0.0233 ~ 0.0233, — - = 0.0016 ± 0.0004. 

Values of ratio m Ul /m U2 and m U2 /m U3 calculated from the values of [i\ satisfying the allowed 
values (Eq. (34)) of neutrino mixing | Vmns I are plotted by sequences of dots in FIG. 2. In 
FIG. 2, sequence of dots in middle is obtained in the case 0i = 0°, 2 = 0°, sequence of dots 
in right-hand side corresponds to the case <p\ = —74.9°, 02 = 0.74°, which are determined 
in quark sector analysis. From this figure, it is shown that the allowed values of ratios 
m Ul /m U2 and m U2 /m U3 are restricted, and for these restricted ratios, the value of |Vmns|i3 is 
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determined as follows, 



771 771 

-^i = 0.44~0.56, — = 0.17 -0.23, | Vmns I is = 0.042 - 0.065, 

TTii/2 m 

for (pi = 0°, (p2 = 0° case, 
" : - 0.92 - 0.98,^ = 0.35 - 0.75, |VWs|is = 0.056 - 0.080, (36) 



for 0i = -74.9°, 2 = 0.74° case. 

Although the case of 0i = —74.9°, 02 = 0.74° is interested in present analysis, but the case 
of 0i = 0.0°, 02 = 0.0° is presented for a comparison. In previous our analysis we 
studied the case of 0i = 0.0°, 02 = 0.0°, and obtained the result, m Ul /m u , 2 = 0.36 — 0.49, 
|Vmns|is = 0.04 - 0.06. For the case of 0i = -74.9°, 2 = 0.74°, it is predicted that the 
neutrino spectroscopy is near the quasi-degenerate(QD) spectrum, and |Vmns|i3 is n °t so 
tiny. 

The allowed values of parameters /j,", M3 are determined from the allowed regions of 
m Ul /m u , 2 , m U2 /m U3 in Eq. (36), as 

2u v2 u u u v 

m U3 « -jg- = (0.047 - 0.053)eV, ^ = 0.045 - 0.050, ^ = -0.014 - -0.019, 
Mi Hq /xg 

^ = ±(0.034 - 0.047), for the case 0i = 0.0°, 2 = 0.0°, 
Mo 

2u u2 u u u u 

m U3 « -t£- = (0.05 - 0.07)eV, ^ = 0.035 - 0.038, ^ = -0.001 - -0.007, (37) 

Ml /Xq /Xq 

^ = ±(0.005 - 0.023), for the case 0i = -74.9°, 2 = 0.74°. 
Mo 

2u u2 u u M 

From the result « 0.06eV, ^ « 0.03 in Eq. (37) and — = 0.0016 in Eq. (35), and 
Mi /Xq M 

assumption /xg - /x ^-f « 73.3GeV in Eq.(19), Mi, M are estimated as 
Mo 

Mi w 1.6 x 10 n GeV, M « 10 14 GeV, (38) 

which is compatible with the result in GUT, Mgut — 2 x 10 16 GeV. For the case 0i = 
—74.9°, 02 = 0.74°, Dirac CP-violation phase 5 is produced and this phase generates CP 
violation in neutrino oscillation. The magnitude of CP violation in v\ — > vy and V\ — > vy 
oscillations is determined by the invariant measure of CP violation, J, as the same definition 
as in Eq. (18'), 

J = Im [(Vmns)i2(VmNs)23(^MNs)i3(^MNs)22] j 
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where (Vmns)i? is a (hj) element of matrix in Eq. (28) without the Majorana phase matrix 
Pm, because Majorana phases in Pm do not contribute in neutrino oscillation [28]. In our 
model, neutrinos are assumed to be the Majorana fermions produced through the see-saw 
mechanism. Majorana nature is found in neutrino-less double beta decay, (/3/3)oi,-decay. The 
(/3/9)oy-decay is characterized by the effective Majorana mass | <m> | defined as 

I <m> | = |m Vl (y M Ns) \l + m V2 (VmNs) 12 + m ^3(^MNs)l3| j (39) 

where, Vmns is a matrix expressed in Eq.(28) with Majorana phase matrix Pm- We estimate 
the Dirac phase 5, the Majorana phases /3,7, the invariant measure J of CP violation, and 
effective Majorana mass | <m> | for the allowed values of A^/Vcm f^/f^o denoted in 

Eq.(37), 

5 = 180.0°, /3 = 0.0°, 7 = 0.0°, J = 0.0, | <m> \ = 0.0059 ~ 0.0079, 

for the case fa = 0.0°, 02 = 0.0°, 
5 = (65.2 ~ 84.3)°, (3 = (24.3 ~ 44.2)°, 7 = (16.9 ~ 31.8)°, (40) 
J = -(0.010 ~ 0.017), I <m> I = 0.026 ~ 0.048, 

for the case (j) X = -74.9°, 2 = 0.74°. 

Thus, for the neutrino mixing derived from the neutrino Dirac mass with (pi = —74.9°, <p2 = 
0.74°, the magnitude of CP violation in neutrino oscillation expressed by J, is predicted as 
to be rather larger than that of quark, which is shown in Eq. (18). Further the effective 
Majorana mass | <m> | is predicted to be not so tiny. 

V. CONCLUSION 

We assumed that the weak bases of flavors (u, c), (d, s), (e, /i), and Dirac neutrino (u e , u^) 
are 53 doublets and t, b, r, and v T are S3 singlets. Further, we assumed that the Higgs 

doublet (i?i,i? 2 ) and Higgs 5*3 singlet Hs- From these 5*3 doublets and singlets, we 
constructed S3 invariant Yukawa interactions and the mass matrices for weak basis of flavors. 



Obtained mass matrices for quark sector are 











\ 




( rt + t% 


A/^e"^ 1 A/^e"^ 1 


\ 


M d = 




/if - /i^e^ 2 








A/^e-^ 1 


li\ - n%e~^ 2 ^e~^ 2 






, A/^e^ 1 


^e^ 2 


l4 


) 




, A/i^e"^ 1 




/ 
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where A = |iif°|/|i?2| and 0i = phase of Hf — phase of H®, and 2 = phase of H® — 
phase of H$. From the present experimental data for quark masses and Vckm matrix in- 
cluding the CP violation phase [271 ] . we can obtain the results Eq. (19), 



Mo 



4.20±0.12GeV, = 0.0120 ± 0.0030, 
Mo 

^ = ±(0.0282 ±0.0008), 
Mo 



— A = -(0.0136 ±0.0004), 
Mo 



Mo 



171.3 ±2.3GeV, 
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Mi 
Mo 



0.00369 ±0.00003, 



M£ 
Mo 1 



-0.00378 ± 0.00003, 



= T (0.0127 ± 0.0007), (opposite sign for that of the ratio /4/yUo) 
Mo 



A = 0.207 ±0.004, 



•(74.9 ±0.8)°, 



(0.74±0.31) c 



CP-phase in standard expression for Vckm is 



5= (68.1 ±3.8)°. 



In our model, the origin of the Cabibbo angle can be explained by the ratio A = \Hi/\H2 \ 
and the origin of the CP violation by the phase difference <fti = phase of H® — phase of Hg. 
For lepton sector, mass matrices are obtained as 



Mr 



M 



D 



A/4^ 1 
y X^e^ 1 

K + M2e" i02 
A/^e"^ 1 
A/^e"^ 1 



A/i^e^ 1 A/4^ 1 \ 



H\ - fi l 2 e i<t>2 



\i\e^ 2 
Mo 



A/i^e" 



■KPl 



A/i^e 



-i<f>l \ 



\ 



Mo 



M, 



M 



( 



Mi 
Mi 
M r 



\ 



(41) 



Diagonalizing the mass matrix Mi and M* f = — M^M^l M^ T obtained through the see- 
saw mechanism, we obtain the charged lepton mass, neutrino mass and neutrino mixing 
Mmns- h- 1 our model, tri-bimaximal-like mixing character of the neutrino mixing Vmns is 
obtained dynamically from the hierarchy of mass parameters fi^ + fi^ ^ Mi ~~ M2 ~ IM3I ^ Moj 
Mi <C M and smallness of A ~ 0.21, without any other symmetry restriction than S3. From 
the present experimental data for charged lepton mass 27| and neutrino mass and mixing 
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Vmns [12|], we obtained the allowed values for mass parameters and |H 



MNS 13; 



A = 0.207 ±0.004, 0! = -(74.9 ±0.8)°, 4> 2 = (0.74 ± 0.31)°, 

(in-put data determined from quark sector analysis,) 

/4 = 1776.84 ± 0.17MeV, = 0.0315, = -0.0324, = -0.0233 ~ 0.0233, 

/i fx /i 

Tfl Tfl 

m ua w (0.05 ~ 0.07)eV, — ^ = 0.92 ~ 0.98, — ^ = 0.35 ~ 0.75, 

IVmnsIis = 0.056 ~ 0.080 ^ = 0.035 - 0.038, ^ = -(0.001 - 0.007), 

^ = ±(0.005 - 0.023), Mi « 1.6 x 10 n GeV, M w 10 14 GeV. 
0o 

Thus, neutrino mass in our model favors a quasi-degenerate (QD) spectrum, and |Vmns| is 
not so tiny. We can estimate the CP violation Dirac phase 8, Majorana phases (/3, 7) and 
CP violation measure J, effective Majorana mass | <m> | in (/3/3)oi/-decay, as 

5 = (65.2 - 84.3)°, p = (24.3 - 44.2)°, 7 = (16.9 - 31.8)°, 
J = -(0.010 ~ 0.017), I <m> I = 0.026 ~ 0.048. 



Appendix: Diagonalization of complex symmetric matrix 

For an arbitrary complex matrix M, M can be diagonalized by the unitary matrix V and 

U as 

V'MU — diag[mi,W2, • • • ,m n ], m,; > 0, (A.l) 

where VtfcfAftV = diagf raf , m\, ■ ■ • , m^}, and we assume that m, 7^ for i ^ k. Since M is 
a complex symmetric matrix, M = M T = £7^ T diag[mi, 777.2, • • • , mjV^, then / 1 M M I ' ' = 
diag[m?, mi, ■ • • , m*]. Hence fdiagK, m|, - • • , m*]V+ = W T dmg[ml m|, • • • , m^][/ T 
This relation implies that 

U T Vdiag{mj, mj, ■ ■ ■ , m^] = diagfm^, mj, ■ ■ ■ , m%\U T V (A. 2) 

Since ?7 T y is unitary matrix, it follows from Eq. (A2) that 

U T V = Sl S = diag[e iQl , e ia \ • • ■ , e ia "}, (A.3) 
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where are arbitrary real constants. Then we can obtain V = U*S^ and SU T MU = 
diag[mi, m 2 , ■ ■ ■ ,m n ]. Thus, we obtain a statement that if M is a complex symmetric 
matrix, M can be diagonalized by the unitary matrix V and U as 

V^MU = diag[mi,m 2 , • • • ,m n ], m { > 0, V^MM^V = diagfrn?, m\, ■ ■ ■ , ml], 
V = U*S\ S = diag[e iQl , e ia \ ■■■ , e ian ], a { : arbitrary real constants. (A.4) 
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